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We investigate the temporal dynamics of Doppler cooling of an initially hot single trapped atom 
in the weak binding regime using a semiclassical approach. We develop an analytical model for the 
simplest case of a single vibrational mode for a harmonic trap, and show how this model allows 
us to estimate the initial energy of the trapped particle by observing the fluorescence rate during 
the cooling process. The experimental implementation of this temperature measurement provides 
a way to measure atom heating rates by observing the temperature rise in the absence of cooling. 
This method is technically relatively simple compared to conventional sideband detection methods, 
and the two methods are in reasonable agreement. We also discuss the effects of RF micromotion, 
relevant for a trapped atomic ion, and the effect of coupling between the vibrational modes on the 
cooling dynamics. 
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Laser cooling of trapped neutral atoms and atomic ions 
is a well established technique: for example, cooling to 
the motional ground state [TJ [21 |S] and motional state 
tomography j4j are routinely performed with resolved 
motional-sideband excitation techniques. Sideband tech- 
niques require the natural linewidth T of the cooling tran- 
sition to be small compared to the vibrational frequency 
of the trapped particle, in order to allow the motional 
sidebands to be resolved. Many experiments are, how- 
ever, conducted in the "weak-binding regime" , where F 
is larger than the oscillation frequency. Here, the cool- 
ing process is essentially the same as Doppler cooling 
of free atoms, because the spontaneous decay process 
is short compared to the atom's oscillation period [S]. 
Even in experimental setups that implement sideband 
techniques, an initial stage of such "Doppler cooling" is 
often employed. The first examinations of Doppler cool- 
ing of trapped ions [3 [i H [M HI did not take 
into account the effects of micromotion due to the trap- 
ping RF field. After cooling and heating effects related to 
micromotion were observed, these effects were explained 
theoretically [TH |TH1 [H] by including the effects of mi- 
cromotion. 

Here, we consider Doppler cooling of a single trapped 
atom or ion. While most previous work has focused on 
the final stages of cooling, our focus will be on the tem- 
poral dynamics of the cooling, particularly in the "hot 
regime" where the Doppler shift due to atom motion is 
comparable to or much larger than F. For the 1-D case 
we find that the cooling rate can be calculated analyti- 
cally in the weak-binding regime without assuming the 
atom to be in the Lamb-Dicke regime. For a trapped ion, 
when we take RF micromotion into consideration, stable. 



highly excited states emerge when only one mode is con- 
sidered [in]. When all three vibrational modes of the ion 
are considered, we find that couplings between the modes 
tend to break the stability of such points allowing cooling 
to reach the Doppler limit. 



A practical application of our results is to estimate 
the initial motional energy of an atom or ion from obser- 
vations of the time dependence of the fiuorescence dur- 
ing the cooling process. As mentioned above, sideband 
spectroscopy is the conventional technique for character- 
izing motional states, and it has been used to charac- 
terize the heating rate of ions in the absence of cooling 
O d [ini [T71 UHl [inj HO] . However, it is more comphcated 
to implement experimentally than Doppler cooling, re- 
quiring more laser beams. Currently, considerable effort 
is being devoted to understanding the anomalous heat- 
ing observed in ion traps [IHl HZl [TSJ [HI H^. A less 
complicated technique for measuring temperature could 
simplify this work. 



This paper is structured as follows: In Sec.|T]we present 
a semiclassical model of the Doppler cooling process for a 
bound atom in the weak-binding regime. In Sees. [ll] and 
|III| we analyze the fiuorescence predicted by the model 
for a single vibrational mode unaffected by micromotion. 
Here, we consider single cooling trajectories and average 
over these with a given distribution of initial motional 
energies. We derive expressions useful for estimating ini- 
tial temperature from fiuorescence observations in these 
sections. Sec. IIVI discusses how to minimize the total 
measurement time required to estimate the mean initial 
energy. In Sec. |V] we consider the effects of other mo- 
tional modes with and without taking into account any 
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RF micromotion experienced by such modes. Sec 
suggests modifications to the basic experimental protocol 
that might provide improved sensitivity of the tempera- 
ture measurements. Sec. |VII| concludes the paper. 
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I. MODEL 



emission, will cause heating at a rate [TT] 



We consider a semiclassical model of Doppler cooling 
of a single weakly trapped atom [11]. We will ini- 
tially consider only a single mode of motion, taken to be 
along the z direction. We assume a harmonic potential 
with oscillation frequency ui^. In Sec. |V]we consider a 
more detailed model that includes three dimensions and 
micromotion for ions. 

The atom is Doppler-cooled by a single laser beam of 
angular frequency J^Lascr and wave- vector fc, detuned by 
from the resonance frequency v„^ of a 



two-level, or "cycling" , transition between two internal 
states, |g) and |e), of the atom. We write the coupling 
Hamiltonian as 



Rabi 



cos(fe 



I'Lasori), (1) 



where x is the atom position, 1i:h is Planck's constant, 
and i^Rabi is the resonant Rabi frequency. 

We assume the atom is weakly bound in the z direc- 
tion, that is, LOz is much less than the excited state decay 
rate F. The atom's level populations are then approxi- 
mately in steady state with respect to the instantaneous 
effective detuning, Aoff = A + A/5 , including the Doppler 
shift Ao = —kzVz, where Vz and kz are the z-components 
of the velocity and wave-vector. The excited state popu- 
lation is then ED 



Pcc{Vz) 



s/2 



l + S+(2Aeff/r)2 



(2) 



Here s is the saturation parameter, proportional to the 
cooling beam intensity, s = 2|riRabiP/r^- 

The excited state population is associated with the 
photon scattering rate dN/dt by the relation dN/dt = 
rpco(wz)- While the momentum kicks associated with 
photon emission are assumed to average to zero over 
many absorption-emission cycles, the absorbed photons 
will impart a velocity dependent momentum transfer due 
to the scattering that can be described by a velocity- 
dependent force 



Fz{Vz) 



dvz 

i— — 
dt 



hkz Tpcc{Vz), 



(3) 



where m is the atom's mass. This velocity-dependent 
force will in general change the motional energy E of the 
atom. If the relative change in energy over a motional 
cycle is small, we can average the effect of Fz over the 
oscillatory motion to find the evolution of E: 



dE 
~dt 



= {VzFz{Vz)) 



(4) 



where the average is over one motional cycle. The average 
energy change per scattering event is dE/dN = HkzVz — 

In addition to Fz, the atom will experience a stochas- 
tic force due to photon recoil that, assuming isotropic 



dE 
~dt 



recoil 



4 [hkz) dN 
3 2m dt ' 



(5) 



where {hkzY /2m is the recoil energy associated with the 
scattering. We will mostly ignore the effects of recoil 
heating in what follows since it will be important only 
near the cooling limit. 



II. ANALYSIS 

We will now analyze the time-dependence of the atom 
fluorescence during the Doppler cooling process, as pre- 
dicted by the model introduced above. 

To simplify the algebra, we will scale energies by h 
times half the power-broadened linewidth, and time by 
the resonant scattering rate: 

{e,5,r} = {E,hA}^}lE,, E,^^Vl + ~s 
2m 2 



T = t/to, 



r 



s/2 



As an example of typical values, we consider a trapped 
25]y[g+ jQjj^ where the — ^ P3/2 cooling transition at 

279.6 nm has a natural hnewidth of T = 27r x 41.4MHz. 
At a detuning of A = -27r x 20 MHz with s = 0.9 and 
kz/k — 0.71, we find that Eo/ks = 1.4mK, where ks is 
Boltzmann's constant, and to = 16ns. The detuning and 
recoil parameters are 5 = —0.70 and r = 0.0018. 

The maximal change in energy per scattering event 
at a given energy is Sm = 2-Jer. The energy at which 
the maximal Doppler shift, which in the scaled units is 
equal to Sm, is equal to the power-broadened linewidth, 
(5j\/ ~ 2, is of interest during the coohng process. For 
reference we note that this energy corresponds to 



1 



£^0 = (1 + s) 



2m 

4 w^: 



(7) 



For the typical experimental parameters considered 
above, E^/r is equal to ks x 700 mK or 3700/11^2 for 
= 27r X 4.0 MHz. 

For harmonic oscillations, the instantaneous Doppler 
shift 8d = ^^d/Fq is distributed according to the prob- 
ability density 



Jo 



fciraclf^D - <^M sin((/))) 



2-K 







if \8d\ < S_ 
otherwise. 



(8) 



where (^oirac is the Dirac 6 function. Since the average en- 
ergy change per scattering event is —6d, and the instan- 
taneous scattering rate is l/{l + 6'^g), where Scs = 5 + 5d, 
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FIG. 2: The scattering rate, dN/dr, vs. time during Doppler 
cooling of a single atom with initial energy e = 10/r for differ- 
ent laser detunings. For 5 < 5c = — 1/\/3, a maximal scatter- 
ing rate occurs at e = Es, as given by Eq. ( |14[ ). The maximal 
value of the scattering rate is given by Eq. (15 1. Closer to 



resonance, the scattering rate increases monotonically during 
the cooling. 



FIG. 1: Scattering rate (top) and energy (bottom) as a func- 
tion of time during the Doppler cooling of a single atom from 
an initial energy of 10/r, as given by Eqs. (lOal and (111 at 
a detuning of half the power-broadened linewidth (5 = —1). 
The dotted curve in the lower plot shows the energy pre- 
dicted by the asymptotic approximation (121. Insets in the 



upper plot show (at two different times) the two components 
appearing in the integral defining dN/dr in Eq. (Ill: The 
probability density of the effective detuning Po {5m ScS — 5) 
(solid) and the Lorentzian line profile, L{5es) = 1/(1-1- 
(dashed) as a function of 5cfi = 5 + Sd at e = 8/r and e — Ss- 
Scattering events where the atom is moving towards the laser 
so that 5d > 0, corresponding to the rightmost peak of Pd, 
result in cooling, and vice versa. The energies of maximal 
cooling and scattering rates, Sc and Ss, are given by Eqs. (131 
and ((iSl). 



the rate of change of e averaged over the secular oscilla- 
tions by Eq. (|4| takes the form: 



de 
dr 



1 



l + id + So^ 



dSr 



(9) 



as illustrated in Fig. [T] We can evaluate the integral as 
detailed in Appendix|A] to find that 



ds 
d^ 



1 



6 

2 Jet' 



(Re(Z) + 5Im(Z)) 
£> (l + (52)/r, 



(10a) 
(10b) 



where Z — Z{6, Sm) ~ W l — (S + i^/Aer. The asymp- 
totic approximation (10b) corresponds to approximat- 



ing Pd{Sm', Sd) by Pd{Sm', 0), which is reasonable in the 
"hot" regime, where the peaks of Pd have small overlap 
with the Lorentz line profile. 

The scattering rate averaged over the motion is analo- 



gous to Eq. ( |10| and is given by 
Pd{5m\ Sd) 



dN 



1 + {S + 6d)'' 



d6r 



1 



2^ 



Im(Z), 



(11) 



as illustrated in Fig. [T] In the limit of £ ^ (1 + S'^)/r, we 
find dN/dr « l/(2-y/er), so that according to Eq. (10b I 
we have in this limit de/dN « S. This corresponds to 
each photon on average extracting an energy of HA. This 
can be understood by noting that in the limit of e » (1 + 
S'^)/r, Pd{Sm]Sd) is to a good approximation uniform 
over the Lorentz line profile, and so the value of Scs — 
5 + Sd averaged over the scattering events will be zero. 
Since each scattering event extracts an energy of —Sd, 
the average cooling per scattering event should indeed 
be S. 

The time-dependence of e is formally found by inte- 
grating de/dr a s giv en by Eq. (10 1. For the asymptotic 
approximation (10b I we find 



e{r) 



-3/2 



e»(l + <52)/r, (12) 



where is the energy at r = 0, as plotted in the lower 
part of Fig. [l] For the exact expression (10a), we must 
resort to numerical methods, although we do find ana- 
lytically that the cooling rate is maximal for e related to 
Shy 



1 



+ S^ ,1 ,1 

cos(- arccos( — 

2r 6 1 



(52 



)), 



(13) 



which quantifies our previous observation that l/r is a 
typical energy scale of the cooling process. 
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The behavior of dN/dr is qualitatively different for 5 
being smaller or larger than a critical detuning, 5c = 
— 1/a/3. For (5 < (5(7, dN/dr has a maximum at 



es= -{6-V3)iS+l/V3). 
4r 



(14) 



For the example parameters listed below Eq. ^ , 5c cor- 
responds to a detuning of ScEq/H — ~2tt x 16.5 MHz. 
Closer to resonance, i.e., when 5c < 5 < 0, no maximum 
occurs, as illustrated in Fig. [2j The maximal scattering 
rate is reached when one of the peaks of the Doppler dis- 
tribution ([s]) is in resonance with the cooling transition, 
as illustrated in the insets of Fig. [T] In the regime where 
a maximum exists, the maximal scattering rate is found 
to exceed the steady state scattering rate by a factor of 



dN 
d7 



dN 
d7 



V3V^l + 5^ 



(15) 



We emphasize that the only approximations made 
above are the weak-binding approximation and the omis- 
sion of recoil heating. In particular, the trapped particle 
is not assumed to be in the Lamb-Dicke regime. For the 
weak-binding regime, 5m > 1 implies that the motion is 
well outside the Lamb-Dicke regime. To find the cooling 
rate predicted by (10 1 in the Lamb-Dicke limit, we note 
that to first order in e we have de/dr w 4fcr/(l + 5'^)'^. 
This corresponds to e decreasing exponentially with r. 
Except for the omission of recoil heating, the value of 
the decay time agrees with previous work that assumed 
the atom was in the Lamb-Dicke regime [51 [52] . 

In the above analysis, we have ignored recoil heating 
as given by Eq. (|5|. In the hmit oi e '3> £c, the ratio of 
heating to cooling is seen to be 4r/(3|(5|), which is a small 
fraction for realistic parameters. For s < the cooling 
is less efficient and the contribution from recoil becomes 
more significant, leading to a nonzero steady-state en- 
ergy. Nevertheless, ignoring recoil heating is reasonable 
when considering only fluorescence, since the scattering 
rate has almost reached its steady-state value when the 
effect of recoil becomes important. We have omitted re- 
coil in this analysis to make 5 the only free parameter 
and simplify the discussion. Recoil can be included in 
calculations by combining Eqs. ([5|, ( |10a[ ), and (111. 



III. THERMAL AVERAGING 

An application of the analysis presented above is to 
estimate the initial motional energy of a trapped atom 
from the fluorescence observed during the cooling pro- 
cess. Using this method, we can estimate the average rate 
of heating experienced by a trapped atom in the absence 
of cooling by first allowing the atom to heat up without 
cooling for a certain period and then observing the time 
dependence of the fluorescence as the atom is re-cooled. 
As discussed in Sec. |ll] we have for e ^ (1 -I- 5'^)/r that 
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FIG. 3: Thermally averaged scattering rate vs. time, for 
5 — — \/3 (top) and 5 = — 1/\/3 = 5c (bottom). In both 
figures, Pq is assumed to be a thermal distribution with e 
equal to 1,2, and 4 times for the solid, dashed, and dash- 
dotted lines respectively. Note that for 5 < 5c, the initial 
fluorescence is larger than the steady-state fluorescence for 
low values of e\ this is attributed to the local maximum in 
the fluorescence illustrated in Fig. [2] 



the average cooling per scattering event is 5. The ap- 
proximate total number of photons scattered during the 
cooling of an atom with initial motional energy Einitiai 
can consequently be approximated by |einitiai/<5|. For the 
example parameters given in Sec. [ll] this corresponds 
to « 800 photons for einitiai = 1/'', corresponding to 
ks X 700 mK. With typical photon detection efficiencies 
of less than 10^'^, very few photons are registered in a 
single experiment. We must therefore repeat many ex- 
perimental cycles consisting of a heating period and a 
cooling period. 

We now consider the form of the fluorescence sig- 
nal when averaged over many such experimental cycles. 
Here, we will assume the heating is stochastic and take 
the distribution Pq{£) of the motional energies at the 
beginning of each cooling period to be the Maxwell- 
Boltzmann distribution with mean energy £, 



Po(e) = |e-^/^ 



(16) 



However, the results below hold for any form of Po(£)- 

The thermally averaged scattering rate is conveniently 
written in terms of the propagator, S, of e: Let 'E.^Eq^t) 
denote the energy at time r of an atom with initial energy 
e(r = 0) = Eq. We can then write the thermally averaged 
scattering rate at time r as 



me'} f 

dr 



de'. 



(17) 



e=S(e',r) 



This can be efficiently computed numerically by noting 
that S(S(£, Ti), T2) = S(£,ri -I- T2), as detailed in Ap- 
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FIG. 4: Experimentally observed fluorescence during Doppler 
cooling of a single Mg^ ion compared with the fluorescence 
predicted by the simple 1-D model. Data points (triangles) 
indicate the observed scattering rates, obtained by integrat- 
ing over many experiments. In each experiment, the time- 
resolved fluorescence is recorded while the ion is re-cooling 
after having been allowed to heat up for a period of 25 s. The 
experimental parameters were those given after Eq. ||6|. Er- 
ror bars are based on counting statistics. The solid curve is 
the scattering rate predicted by Eq. (17 1, assuming the mo- 



tional energy of the ion after the heating period to be given 
the Maxwell-Boltzmann distribution ( 16 1 with e = 5.1/r, cor- 
responding to a te mp erature of 3.9 K. Since e is the only free 
parameter of Eq. (171, the estimated value was extracted by 



a single parameter fit, and agrees reasonably well with an in- 
dependent temperature estimate of 3.4 ± 0.3 K extrapolated 
from heating rates measured in the same trap by use of the 
Raman sideband technique [20j . 



pendix |B] Figure [3] shows the thermally averaged scat- 
tering rate for a few different parameters. 



The fluorescence predicted by Eq. (17) has been found 
to be in good agreement with experimentally observed 
fluorescence. We show one experimental data set for 
comparison in Fig. |4] the experiments are more fully 
described in [5D]. Furthermore, the resulting estimated 
heating rates have been found to agree well with results 
obtained using the Raman sideband technique [ITl |2Q] . 
This agreement may at first seem surprising, given that 
the two methods probe very different energy scales. For 
the measurements based on the Raman sideband tech- 
nique the ion was only allowed to heat for a few mil- 
liseconds thereby gaining a few motional quanta while 
the measurement results presented in Fig. |4] are based on 
25 s heating periods allowing the ion to gain many mo- 
tional quanta. However, the results should agree if, as 
expected, the heating rate is constant over these energy 
scales. 



IV. OPTIMAL EXPERIMENTAL 
PARAMETERS 

We now examine how the total measurement time re- 
quired to reach a given accuracy on the heating rate esti- 



FIG. 5: Total measurement time required to establish the 
heating rate with a given accuracy, assuming the integration 
time to be dominated by time for re-heating. A local maxi- 
mum is observed for all detunings, but for & — — 1/-\/3 it is 
located outside the range covered by the plot. 



mate depends on the choice of experimental parameters. 
As the recoil parameter r will be fixed by choice of atom, 
we consider only the choice of optimal values for e, b, and 
laser beam intensity. 

It is clear from Fig.[2]that the relevant size of the signal, 
in terms of fluorescence photons emitted, for a given ini- 
tial motional energy increases with decreasing detuning: 
the re-cooling is slower and the change in scattering rate 
is larger. For a given experimental setup, the optimal 
detuning is decided as a compromise between re-cooling 
signal and ability to re-cool atoms that have been highly 
excited by e.g., collisions. 

For a given value of e, the experimental signal, in terms 
of the number of photons scattered before steady state 
is reached, does not depend on the laser beam inten- 
sity. Since e is the average initial energy relative to Eq, 
which is proportional to the power-broadened linewidth, 
a lower laser beam intensity will give a larger signal for 
a given heating period. This suggests using the smallest 
feasible laser intensity, requiring a compromise with re- 
spect to robust cooling and detector dark counts. From 
this standpoint we want to keep the saturation parameter 
below, but probably close to, 1. 

For a given detuning and laser intensity, an additional 
choice of the length of the heating period in each ex- 
perimental cycle has to be made: Should we perform a 
relatively low number of cycles with long heating periods 
or more cycles with shorter heating periods? 

To answer this question, we estimate the total mea- 
surement time, Ttot, required to reach a certain relative 
accuracy on the estimate of the heating rate. We assume 
a constant heating rate and assume that the total time is 
dominated by the heating periods, so that Ttot is propor- 
tional to the average initial energy, E oc £ \/l + s, and to 
the the number of runs. 

We will consider a setup where the observed fluores- 
cence is collected in sequential time-bins that are short 
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compared to the total time required for the cooling pro- 
cess. In the limit where the distribution of the integrated 
number of counts, rii, in time-bin i is described by a nor- 
mal distribution with variance ai, we can estimate the 
uncertainty on the maximum-likelihood estimate of e for 
a given dataset by |23] 



~di 



(18) 



It follows from Eqs. (lOal, (111, and (17l, that in the 1 



D case the cooling dynamics can be rewritten in a form 
independent of r by reparametrizing in terms of A^r, er, 
and TT. We will denote the reparametrized scattering 
rate by 



(19) 



Since the relative uncertainty on the heating rate esti- 
mate is equal to a(e) je and ui — ^/ui, Eqs. ( 18 ) and ( 19 ) 



allow us to estimate the time required to obtain a given 
relative uncertainty on the heating rate: 



Ttot 



er 



/ dRs{er,q) 
\ der 




(20) 
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FIG. 6: The probability density (221 for the combined 
Doppler shift 5'^ = 



due to two excited modes 



(curve on left wall) is a marginal distribution of the joint prob- 



ability density Po{5^^\5^^'' 



)Pd{5 



(v) Av) 

AI ' "d 



) Of i5^^\5^^^) (3-D 



surface) . Note that only two of the peaks in the joint proba- 
bility distribution lead to peaks in the marginal distribution. 



Plot is drawn for S^^j' = 3 and S)^/ = 1, as the dash-dotted 



(y) 



\5^jj\ < 0.95 5)^.1 for illustrational purposes. 



line in Fig. 7(a) and the joint distribution is truncated to 



Note that the right hand side depends only on S and er. 
Fig. [5] shows Ttot + s calculated for different detun- 
ings. The figure confirms that a low detuning is indeed 
favorable, and also shows that for a given detuning, Ttot 
decreases with increasing e. This is not surprising, given 
that the time to cool by a certain amount of energy in- 
creases with atom temperature, as illustrated by Fig. [T] 
It is clear from Fig. [sjthat the heating period should be 
chosen long enough to get a significant signal, i.e., e > Ec, 
but the optimal heating period must be decided based on 
other experimental parameters such as trap depth and 
background gas collision rate. 



V. COOLING IN THREE DIMENSIONS 

So far, we have considered only cooling in one dimen- 
sion. In this section we will consider the effect of the 
vibrational modes in other directions on the cooling pro- 
cess. Our goal is to gain a qualitative understanding of 
the effects of the transverse modes on the cooling dy- 
namics of the z mode, with the intent of establishing to 
what extent the simple 1-D model presented above is a 
reasonable approximation. 



A. 3-D cooling of neutral atoms 

For a neutral atom, the confinement transverse to z is 
not associated with micromotion, as it is for ions, and the 



1-D weak-binding model extends immediately to three di- 
mensions. Let Ei, i = {x,y,z}, denote the motional en- 



ergy in mode i, 5 



(i) 



-hkiVi/Eo the Doppler shift, and 



6\j = 2y/eiri the maximum Doppler shift. Although all 
modes are formally identical in the absence of micromo- 
tion, we will discuss the cooling dynamics with a focus 
on the z mode. 

In experiments it is typically easy to make the frequen- 
cies of the three modes incommensurate, which we will 
assume here. In that case, we can write the rate of change 
of e, as 



dez 



5'j^^)dS 



(^) 

D ' 



(21) 



where is the effective line profile experienced by the 
z mode, obtained by convolving the Lorentz line profile 
with the distribution P^^'' of the combined Doppler shift 

"d — "d 



due to the x and y "spectator" modes, 



Ph'\s'^) ^ / PD{S^^i^;u)PD{S'i'j>;S'^ - u)du. (22) 



verges) at 5'j-, = ±\5^^ 



H 

('if 

M 



As illustrated in Figs. |6l and 7(a) P^'^ is peaked (di 



If |<5 
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-'M I 



> 2, 



the peaks are separated by more than the width of the 
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FIG. 7: Cooling of the z mode with two excited spectator modes in the absence of micromotion. Part (a) shows the probabiHty 
density P^^' of the combined Doppler shift due to the x and y modes for equal to (0,4) (solid), (2,2) (dashed), and 

(3, 1) (dash-dotted). Part (b) shows the effective line profile /^^((J^g ) obtained by convolving P''-^^ with the Lorentz line profile. 
Here 5^^' = (5-|-(5^^ The spectator mode parameters are the same as for (a), and the dotted line show the Lorentz line profile, 
corresponding to the spectator modes being cold, essentially the 1-D cooling case. Part (c) shows the z mode cooling rate as a 
function of for the effective line profiles of (b). 



Lorentz profile, and Rz will be double-peaked, as illus- 
trated in Fig. |7(b)[ It follows from Eq. (|2l]) that the 



cooling rate in the limit of small b'^^^ is proportional to 
the slope of Rz at 5, and that the rate of change of is 
positive if the slope is negative. If — \b'^^} — < (5 < 0, 
this will result in heating of the z mode, at least as long 

i.e., while 
(5, as illustrated in Fig. 7(c) 



as 5 ± J^^"* are both inside the peaks of R 



for 



(•^M'^m) = (0-4) and (1,3) [HI [E]. The figure also 
shows that this thermalization or energy equilibration 



effect is not present if \5 



5, 1, as Rz is not 



double-peaked in this case. Mathematically, dez/dt, Rz, 

(z) 

and Pj^ are all conveniently expressed as convolution 
integrals of functions with known Fourier transforms. 

The dashed lines in Fig. [8] show the cooling rates pre- 
dicted by Eq. (21 1 for the case of only one excited spec- 



tator mode at different energies. When (5^-* is large com- 

pared to S)^ , we see that the cooling rate is almost un- 
affected by the spectator mode. This can be understood 
by noting that in the limit where Pd (S]^^ ; 6^j^^ ) is uniform 

( z) ( z) 

over the values of S)j where Rz{S + S}^ ) is nonzero, the 
symmetry of Rz implies that the average energy change 
per scattering event isS, as also discussed in Sec. [IT] Since 



Rz{Scs) « for 6^s > 1 + S 



M 



r(y) 



this implies that 



the temperature of the spectator modes will not affect 
the cooling rate in this limit. At lower values of ^'■^^ 



M 



we 



generally see a decrease in the cooling rate in a gradual 
approach to the thermalization regime discussed above. 

The consequences of thermalization/equilibration pro- 
cess are complex, when considering the full 3-D cooling 
problem. Consider for instance the case where only one 
mode is initially hot. According to the discussion above, 
this will result in heating of the two remaining modes. 



until the fastest heating mode has reached a value of Sm 
similar to that of the initially hot mode. After this ther- 
malization, the modes will be cooled simultaneously at a 
cooling rate significantly lower than the cooling rate for 
a single hot mode. 

At this point, it is worth reconsidering the validity of 
our omission of recoil heating: The recoil heating rate 
as given by Eq. ([s]) is seen to have a maximum value 
of 47V3 at the resonant scattering rate. It is clear from 
Fig. [sjthat for typical values of r on the order of lO^'^, 
recoil is insignificant at high energies. 



B. 3-D cooling of ions including the effects of 
micromotion 



For an ion in a linear Paul trap, if we take the z direc- 
tion to be the axis, confinement in the transverse x and 
y directions is provided by the ponderomotive potential 
of an RF quadrupole field. The full 3-D cooling problem 
including micromotion on the transverse modes is very 
complex even in the Lamb-Dicke regime [13l [HI [21] . At 
low saturation, the effects of the micromotion caused by 
an RF field of frequency f2 can be modeled by including 
micromotion sidebands in the line profile [T3]: Micromo- 
tion with peak amplitude a(x(t)), where x is the ion 
position averaged over one period of the RF field, can be 
described by the line profile, 



1 



+ (Jeff - nfi)2 ' 



(23) 



where SeB 



J2i^D^ is the effective detuning, /3 



\a{x) ■ k\ is the micromotion modulation index, f2 
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FIG. 8: Predicted cooling rates for (a) a transverse mode in the presence of an excited z mode and (b) the z mode in the 
presence of a single excited transv erse mode at 5 = —1. Dashed lines assume that the trapping potential has no associated 
micromotion, as described by Eq. (21 1. In this case, the cooling rates are of course identical for the two cases. The solid lines 
assume the transverse confinement to be provided by an RF quadrupole potential with a frequency (7 = 6, as described by 
Eq. ( 24 1 . For both situations, the cooling rate is plotted with the spectator mode having a motional energy of er = {0, 1 , and 4} . 



The plots are based on the weak binding model, and thus assume secular frequencies to be small compared to the linewidth. 
Note that the thermalization/equilibration effect discussed in the text is clearly observed for both cases. 



hKljE^ is the scaled RF frequency, and J„ is the n-th 
Bessel function. 

In contrast to the situation in Ref. [13], we are con- 
sidering a case where (5 changes during the secular mo- 
tion. Since /3 and Jeff depend on x and S, respec- 
tively, we parametrize the secular motion by the in- 
stantaneous phases, 0^, where Xi = cos(0i(t)), and 



with a micromotion sideband, as described by Eq. ( 23 1 



where 



,(0) 



is slowly varying and 



Choosing 



the X and y axes so that the RF field is proportional 
to (xx — yy)cos(rit), we find that (5^^ — s'"^] sm{(f>i). 
In the limit where the transverse confinement is modi- 
fied only weakly by static potentials, so that uj^ « LUy, 
we find in the pseudopotential approximation that /3 = 



cos((/)3;) — S^j cos{(j)y) \ /Cl, which we note to be 
independent of the secular frequencies. In this case we 
have 



dT 



- J S,R^{S,s{cl,),P{<t>..<f>y))-^, (24) 



where the integral is over [0, 27r] in all dimensions. Note 
that since the modulation index depends only on the 
transverse components of the motion, the effect of ex- 
cited transverse modes on the cooling of the z mode can 
still be described in terms of an effective line profile, as 
in Eq. ([2ll. 



For the cooling of the transverse modes, the effects of 
micromotion on the cooling rates is pronounced, as il- 
lustrated by Fig. 8(a) A very clear qualitative difference 



from the cooling rate in the micromotion-free case is that 
at sufficiently high RF frequencies {(l > 4.4 for S — —1), 
stable points for the transverse mode energies develop 
even when the remaining modes are cold. This effect has 
been discussed in Ref. flS^ and is attributed to the heat- 
ing peak of the Doppler distribution becoming resonant 



This might be related to the bistable behavior reported 
in some single ion experiments [24l |25l [26] . The stability 
breaks down when thermalization is taken into consider- 
ation. Consider for instance the stable point indicated 
by Fig. 8(a) to exist for e « (IS/r^;, 0, 0). Here, it is 
clear from the figure that when the z mode has heated 
to Sz > I/tz, cooling of the x mode will commence. 

When 5m ^ \/2 for the transverse modes, we find 
that only the Jq term of Eq . ([23| contributes significantly, 
and the argument of Sec. |V A| that the cooling rate for 
the z mode is not affected by excited transverse modes 



when jI^'' > 1 



also applies here, as illustrated 

by Fig.'[8(b)] 

It is clear from the results above that we cannot ignore 
the transverse modes if their associated maximal Doppler 
shifts are comparable to that of the z mode. If, how- 
ever, we assume the transverse modes are cold enough to 
avoid the heating effects described in Figs. [7] and [8] we 
have seen above that the primary effect of the transverse 
modes will be to slow down the cooling of the z mode. 
This would result in the 1-D model overestimating the 
mean initial energy of the z mode. However, for many 
experiments that use linear RF traps, it is reasonable to 
assume that the transverse modes are heated significantly 
less than the z mode. This is because most investiga- 
tions of the anomalous heating in ion traps have found 
the results to be consistent with heating rates having a 
frequency dependence of w"" with n > 1 [16l [HI [20] . 
Since the transverse mode frequencies are often an order 
of magnitude larger than w^, this would indeed lead to 
the transverse modes being significantly colder than the 
z mode. Also, since the energy in the transverse modes 
only affects the cooling of the z mode through the re- 
sulting Doppler shift, the effect of the transverse modes 
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FIG. 9: Averaged fluorescence vs. time as predicted by 
Monte Carlo simulations for best fits to the data presented 
in Fig. |4] for different models of the frequency dependence of 
the anomalous heating rate: The solid line assumes transverse 
modes to be unaffected by the heating; this should be very 
well approximated by the best fit to the 1-D model (dotted), 
as also plotted in Fig. [4] For the dashed and dash-dotted 
lines, we assume heating to be proportional to oj"^ '* and in- 
dependent of uj, respectively. The fitted value of the z mode 
temperature for the three cases is 3.9 K, 3.7K, and 1.4K 



could be further reduced by aligning k to have a smaller 
projection on the transverse modes. This would however 
reduce the efficiency of cooling of the transverse modes 

Finally, another effect with respect to micromotion is 
that the presence of uncontrolled static stray fields can 
result in the ion experiencing micromotion even at the ion 
equilibrium position. At temperatures where 5m 
the first order effect according to Eq. ( 23 1 of this will be 



a reduction of the central spectral component by a fac- 
tor of Jo(/3)^; see for example Ref. [2.7. • We note that 
this effect can be compensated by using an effective sat- 
uration parameter based on the steady-state fluorescence 
observed in the trap. 



optical Bloch equations through a large number of cool- 
ing trajectories. For each trajectory, we propagate the 
density matrix p of the ion's internal state according to 
the master equation 



dp 
lit 



,H^^\x,t)\+2LpL^ - {L^L,p] , (25) 



where L = |g) (e| \/Y /2 is the Lindblad operator for 
excited state decay and x(t) — x{t) + a{x{t)) cos{^lt) 
for the X and a introduced above. Coupling to the 
motional state is modeled by the average light force, 
mx = hkT Pcc(t). This model assumes neither the atoms 
to be weakly bound nor the cooling beam intensity to be 
low but does neglect recoil heating. 

Figure[9]shows the result of fitting simulations with dif- 
ferent assumptions for the frequency dependence of the 
heating to the dataset presented in Fig. [4j We find that 
if we assume the transverse modes are not heated, we ob- 
tain a temperature estimate of 3.9 K, in agreement with 
the result of fitting the 1-D model to the data, as il- 
lustrated by Fig. [4] The z mode temperature of 3.7 K 
estimated from the uj^^-^ model is close to, and slightly 
smaller than, this value, and agrees with the temperature 
estimate of 3.4 ± 0.3 K based on extrapolating heating 
rates measured with the Raman sideband technique for 
the same trap configuration. This particular form of the 
frequency dependence of the heating rate was observed 
for the same trap when the Raman sideband technique 
|20j was used, and similar frequency dependencies have 
been observed in other geometries [HI [18] . If we instead 



assume an uj~ dependence of the heating, the results 
only change slightly. 

Our main conclusions from the simulation results are 
that the primary effect of the presence of weakly heated 
spectator modes will be to slow down cooling due to ther- 
malization. If uj~^''^ heating of the transverse modes is 
assumed, the 1-D model will somewhat overestimate the 
motional temperature of the axial mode. 



C. Departures from the weak-binding, 
low-saturation limit 

In most experimental situations, we will not strictly 
fulfill the requirements of low saturation or weak binding. 
In particular, for the trap referenced in Fig.[4]the secular 
frequencies of the transverse modes are approximately 
equal to half the 41.4 MHz linewidth of the Doppler cool- 
ing transition, making the weak binding assumption only 
approximate. Also, the illustrated data were obtained at 
a saturation parameter of 0.9, outside the validity region 
of the line-profile model that accounts for RF micromo- 
tion ( [24| . To validate our claim that the fiuorescence sig- 
nal predicted by the 1-D model is a good approximation 
if the heating rate is assumed to be a strongly decreas- 
ing function of w, we performed a numerical Monte Carlo 
simulation of the fiuorescence, based on integrating the 



VI. MODIFIED EXPERIMENTAL PROTOCOLS 

We consider two modifications to the experimental pro- 
tocol to reduce the total measurement time. Both are 
motivated by the fact that the size of the signal from a 
given amount of heating increases with increased initial 
energy. 

One approach would be to coherently add a known 
amount of energy to the z mode at the start of the heating 
period. If the added energy is enough to bring the atom 
into the slow-cooling regime, this will increase the signal 
change due to a given amount of additional heating, as 
illustrated in Fig. [lO) 

Alternatively, parametric amplification [28j [29l [30] 
could be employed after the heating cycle to modify the 
thermal distribution. Parametric amplification can be 
implemented by modulating the z trap potential at 2wz, 
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FIG. 10: Improving the sensitivity of the temperature mea- 
surement by deliberate excitation. The dashed line shows 



the average scattering rate (171 as a function of time for an 



atom that was deliberately excited to a motional energy of 
So ~ Ssc, so that Po(e) = foirac(£ — eo)- The solid line shows 
the scattering rate as a function of time for an atom which 
has first been deliberately excited to a motional energy of 
£o ~ 5£c and then allowed to heat for a duration which added 
an average thermal energy of e = 0.25/r. For comparison, 
the dash-dotted line shows the signal for an atom experienc- 
ing the same heating period without any initial excitation, 
i.e. with s = 0.25/r. In all cases, S = —1. 
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FIG. 11: Improving the sensitivity of the temperature mea- 
surement by parametric amplification. The dashed and dash- 
dotted lines show the average scattering rate vs. time during 
the re-cooling of an atom which has been allowed to heat to 
£ = 0.9£c and I.Iec respectively at 5 = —1. The solid and 
dotted lines show the average scattering rate for the same 
conditions, except that the motional amplitude has been para- 
metrically amplified by a factor of 3 after the heating period. 



and leads to amplification of one quadrature of the mo- 
tion while damping the other quadrature. For a low value 
of e, parametric amplification would increase the fraction 
of experiments in which the atom is in the slow-cooling 
regime at the beginning of the cooling process, thus in- 
creasing the signal for a given heating period, as illus- 
trated in Fig. 11 



VII. CONCLUSION 

In conclusion, we have shown that the motional en- 
ergy of a trapped atom or ion can be estimated from the 
temporal changes in fluorescence observed when Doppler 
cooling is applied. Specifically, the initial energy can be 
estimated by fitting Eq. (17), where the mean initial mo- 



tional energy is the only free variable, to the observed 
fluorescence. Our analysis assumes the oscillation fre- 
quency of the atoms is much smaller than the linewidth 
of the optical transition used for Doppler cooling and the 
motional energy at the start of the cooling is thermal. 

Compared to Raman sideband transition methods for 
heating rate measurements, this method is simpler to im- 
plement experimentally but requires longer measurement 
duration for traps with low heating rates. On the other 
hand, for high heating rates, where sideband cooling is 
inefficient, this may be the method of choice. We have 
shown that in the typical situation, where the time for 
heating dominates, the total measurement time decreases 
with decreasing laser intensity, decreasing laser detun- 
ing, and increased heating period duration. We have 
compared the trade-off between these par ameters (Fig. [s] 
and Eq. (20)). Finally, we show (Sec. IV I that the total 
measurement time can be reduced by adding additional 
energy to more quickly bring the ion into the low fluo- 
rescence regime. 

By comparison with various models of 3-dimensional 
Doppler cooling, we have established that under typical 
experimental conditions the effects of the high-frequency 
modes are small, and that they will lead to temperature 
estimates that are somewhat higher than the actual tem- 
perature of the low- frequency mode. 
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APPENDIX A: INTEGRALS 

The integrals appearing in Eqs. (|9]), (111, and (21) 
are all convolution integrals of elements with analytical 
Fourier transforms and can thus be easily evaluated in 
Fourier space. For the 1-D integrals, the inverse Fourier 
transform can also be performed analytically. Here we 
present a more direct approach to evaluating the 1-D in- 
tegrals. 

For a, 6 e M we define Z{a, b) as 



Z . 



1 



sin(0) — z 27r 



(Al) 
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where z = {a + i) /b. Noting that 

^ =b ^ ,^ ... iibx~a)+i), (A2) 



X — z 1 + (a — bxy 



we find that, according to Eq. (All 



^'^ 6sin(0) d(j> _ 1 
l + (a-fesin(0))2 2^ ~ 6 

27r 



1 



1 + (a - fesin((/)))2 27r 



(Rc(Z) + aIm(Z)) 



Ini(Z). 



Taking the branch cut discontinuity for \A to be along the 
negative real axis, we have for b > that y/ {—iz)'^ = —iz, 
so that 



Z{a,b) = 



b>0. 



(A3) 



APPENDIX B: NUMERICAL CALCULATION OF 
THE AVERAGED SCATTERING RATE 

In this section we present an efficient numerical method 
for evaluating the averaged scattering rate given by 



Eq. (17) 



Introducing £„ = S(£o, nAr), for n = 0, 1, . . ., we note 
that S(e„j,nAr) = Sm+n- The values of e„ are the en- 
ergies along a single cooling trajectory. If the scattering 
rate can be considered constant on time scales of Ar, 



'^^ (-( W T? 



T e [(n- l)AT,nAr], (Bl) 



we find that the thermally averaged scattering rate, as 
given by Eq. (17), averaged over the same intervals can 



be approximated by 

oo 



m=0 



(B2) 



Since the values of the i?„ are independent of Pqi Rn is 
easily calculated for different Pq by list convolution. 
For the Maxwell-Boltzmann distribution, a numeri- 



cally stable form of the weight factors appearing in (B2| 
is 



6+Ae/2 



e-Ae/2 



-s je 



2e~''/^sinh 



Ae 
2^ 
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